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Abstract 

We perform a global analysis of neutrino oscillation data in the framework of three neutrinos, 
including the recent improved measurement of the neutral current events at SNO. In addition to 
all current solar neutrino data we take into account the reactor neutrino data from KamLAND 
and CHOOZ, the atmospheric neutrino data from Super-Kamiokande and MACRO, as well as 
the first spectral data from the K2K long baseline accelerator experiment. The up-to-date best fit 
values and allowed ranges of the three-flavour oscillation parameters are determined from these 
data. Furthermore, we discuss in detail the status of the small parameters a = Am^QL/Am^^^j and 
sin^ 013, which fix the possible strength of CP violating effects in neutrino oscillations. 
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I. INTRODUCTION 



Recently, the Sudbury Neutrino Observatory (SNO) experiment jl| has released an im- 
proved measurement with enhanced neutral current sensitivity due to neutron capture on 
salt, which has been added to the heavy water in the SNO detector. This adds precious 
information to the large amount of data on neutrino oscillations published in the last few 
years. Thanks to this growing body of data a rather clear picture of the neutrino sector 
is starting to emerge. In particular, the results of the KamLAND reactor experiment 
have played an important role in confirming that the disappearance of solar electron neu- 
trinos the long-standing solar neutrino problem, is 
mainly due to oscillations and not to other types of neutrino conversions 14 , ^ . Moreover, 
KamLAND has pinned down that the oscillation solution to the solar neutrino problem is 
the large mixing ang le MSW solution (LMA-MSW) E3, Q, Q characterized by the 



presence of matter effects j20|j|2l|. On the other hand, experiments with atmospheric neu- 

show strong evidence in favour of z/ ,, — » oscillations, in 



trinos 



231 |24, 



agreement with the first data from the K2K accelerator experiment [28[. Together with the 
non-observation of oscillations in reactor experiments at a baseline of about 1 km [2^ [s^l , 
and the strong rejection against transitions involving sterile states 3l|, these positive evi- 
dences can be very naturally accounted for within a three-neutrino framework. The large 
and nearly maximal mixing angles indicated by the solar and atmospheric neutrino data 
samples, respectively, come as a surprise for particle physics, as it contrasts with the small 
angles characterizing the quark sector. 

In this work we report on the implications of the new SNO-salt measurement on the 
determination of three-neutrino oscillation parameters. In Sec. [H] we investigate the impact 
of the new data on the solar neutrino parameters by performing a fit to solar data and the 
KamLAND reactor experiment in a simple two-neutrino framework. In Sec. IIIII we present 
the results of the general three-neutrino fit to the global data from solar, atmospheric, re- 
actor and accelerator neutrino experiments. After fixing our notation in Subsec. IIII Al we 
investigate the impact of three-flavour effects on solar neutrino parameters (Subsec. Illl B|) 
and we give the current best fit points and allowed ranges for all three-flavour oscillation pa- 
rameters (Subsec. III! C|) . Furthermore, in Subsec. Qll Dl we discuss in detail the constraints on 
the small parameters sin^ ^13 and the ratio ArngQ^/^^^iTM) which are the crucial parameters 
governing the possible strength of CP violating effects in neutrino oscillations. A summary 
of the constraints on three-neutrino oscillation parameters from all current experiments is 
given in Sec. IIVI 
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II. TWO-NEUTRINO ANALYSIS OF SOLAR AND KAMLAND DATA 



In this section we investigate the impact of the new SNO-salt data on the determination 
of the solar neutrino parameters sin^ 6'sol and Aml^^^ in a two-neutrino framework. Before 
presenting the global analysis of all solar data it is instructive to determine the restrictions on 
solar neutrino oscillation parameters that follow from SNO data only. The regions delimited 
by the lines in Fig. ^ are obtained only from the new SNO-salt data in the form of the 
neutral current (NC), charged current (CC) and elastic scattering (ES) fluxes as reported 
in Ref. To include these data we follow the prescription given in Ref. (s^, taking 
carefully into account of correlations and systematic effects. The colored regions of Fig. ^ 
result from the total SNO data including the new salt measurements 0], as well as the 2002 
spectral day/night data [l^. Following Ref. 0] the 2002 and 2003 data are treated 
uncorrelated. For details of our SNO spectral analysis see Ref. j3l|. One sees that Fig. [T] 
is in good agreement with the result obtained in Ref. P|. This figure shows also that other 
experiments play an important role in ruling out low-mass solutions, such as the LOW 
solution, the small mixing angle non-adiabatic branch of the MSW plot, as well as the 
vacuum type solutions, still present here at the 99% C.L.. 
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With this cahbration at hand, we now turn to the full two-neutrino data analysis of the 
erlobal set of solar neutrino data. In addition to the SNO data we take into account the 
rates of the chlorine experiment at tire Homestake mine S 3 ("S ± 0.16 ± 0.16 SNU), the 
most up-to-date results of the gallium experiments SAGE 0, 0] (69.1 SNU) 
and GALLEX/GNO QHll (69.3 ± 4.1 ± 3.6), as well as the 1496-day Super-Kamiokande 
data sample [lO| in the form of 44 bins (8 energy bins, 6 of which are further divided into 
7 zenith angle bins). The analysis methods used here are similar to the ones described in 
Refs. j^, 34| and references therein, with the exception that in the current work we use the 
so-called pull-approach for the calculation. As described in Ref. ji^, each systematic 
uncertainty is included by introducing a new parameter in the fit and adding a penalty 
function to the x^- However, the method described in Ref. j^l is extended in two respects. 
First, it is generalized to the case of correlated statistical errors J36|, as necessary to treat 
the SNO-salt data. Second, in the calculation of the total we use the exact relation 
between the theoretical predictions and the pulls associated to the solar neutrino fluxes, 
rather than keeping only the terms up to first order. This is particularly relevant for the 
case of the solar flux, which is constrained by the new SNO data with an accuracy better 
than the prediction of the Standard Solar Model (SSM) js?]. In our approach it is possible 
to take into account on the same footing both the SSM boron flux prediction and the SNO 
NC measurement, without pre-selecting one particular value. In this way the fit itself can 
choose the best compromise between the SNO NC data and the SSM value. 

In Fig. 121 we compare the allowed regions for the oscillation parameters before and after 
the new SNO-salt data. One finds that especially the upper part of the LMA-MSW region 
and large mixing angles are strongly constrained by the new data. This follows mainly from 
the rather small measured value of the CC/NC ratio of 0.306 ± 0.026 ± 0.024 0, since 
this observable increases when moving to larger values of ArngQ^ and/or sin^6'soL (see, e.g., 
Ref. jjisl). The best fit values for the parameters are 



sin^ 



0.29, 



Aml^^ = 6.0 X 10~^ eV^ (solar data). 



(1) 



The new data show a much stronger rejection against maximal solar mixing: from the Ax^ 
projected onto the sin^ 6'sol axis shown in Fig. |21 we find that sin^ ^sm, = 0.5 is excluded at 
more than 5a, ruling out all bi-maximal models of neutrino mass j38|. 

We now turn to the combination of solar data with data from the KamLAND reactor 
experiment . To this aim we use the event-by-event likelihood analysis for the KamLAND 
data as described in Ref. 



binned data. In Fig. El we show the results o 



39[, which gives stronger constraints than a x^^fit based on energy 
the combined analysis. Comparing post and 



pre-SNO-salt results (see, e.g., Refs. |ii,E3,EiEi) one finds that the new data disfavour 
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Figure 2: Projections of the allowed regions from all solar neutrino data at 90%, 95%, 99%, and 3a" 
C.L. for 2 d.o.f. onto the plane of sin^ ^sql and Arn^Q^ before (lines) and after (colored regions) the 
inclusion of the SNO-salt data. Also shown is Ax^ as a function of sin^ ^sql and Arn-gQ^j minimized 
with respect to the un-displayed parameter. 

the high Am^ region, which appears only at the 99.5% C.L. (A^^ = 10.7).^ The best fit 
point of the global analysis occurs at 

sin^ ^soL = 0.30 , Am^oL = 6.9 x 10"^ eV^ (solar+KaniLAND data). (2) 

We note that for the first time it is possible to obtain meaningful bounds on solar neutrino 
parameters at the 5a level, showing that neutrino physics enters the high precision age. 
From the projections of the onto the AmgQL and sin^ 6'sol axes also shown in Fig. El we 
find the following allowed ranges at 3a (5cr) for 1 d.o.f.: 

0.23 (0.17) < sin^^soL < 0.39 (0.47) , 

(3) 

5.4 (2.1) X 10"^ eV^ < Am^^^ < 9.4 (28) x 10~^ eV^ . 



Recently, an improved determination of the day/night asymmetry in the 1496-day Super-Kamiokande 
data has been released 0|, which leads to a further disfavouring of the high-LMA region. Currently 
not enough information is available to reproduce this result outside the Super-Kamiokande collaboration. 
Therefore, tlie day/night data is treated as in our previous analyses [sif . 
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Figure 3: Projections of the allowed regions from all solar neutrino and KamLAND data at yu7o, 
95%, 99%, and 3a C.L. for 2 d.o.f. onto the plane of sin^ 9sol and Am^Q^ before (lines) and after 
(colored regions) the inclusion of the SNO-salt data. Also shown is Ax^ as a function of sm "sol 
and Attt-Iql) minimized with respect to the un-displayed parameter. 

III. GLOBAL THREE-NEUTRINO ANALYSIS 
A. Notations 



In this section the three-neutrino oscillation parameters are determined from a global 
analysis of the most recent neutrino oscillation data. For earlier three-neutrino analyses 
see Refs. jisl. Is^ 41 1. To fix the notation, we define the neutrino mass-squared differences 
AmgQL = Am|i_= ml — m\ and Am^^,^ = ^rn\i = ml — mf, and use the standard parame- 
terization l^A k^l for the leptonic mixing matrix: 



U 



( C13C12 ■S12C13 Si3 \ 

-S12C23 - S23S13C12 C23C12 - S23S13S12 S23C13 
\ ■S23S12 — S13C23C12 — S23C12 — ■S13S12C23 C23C13 J 



(4) 



where = cos% and Sij = •mvQiy Furthermore, we use the notations Q\2 = ^sol and 6*23 = 
6'atm- Because of the hierarchy AmgQ^ ^ ^'^atm is a good approximation to set AmgQ^ = 
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Figure 4: Sections of the three-dimensional allowed regions in the (sin^ ^sql) '^"'-sol) plane at 
90%, 95%, 99% and 3a C.L. for 3 d.o.f. for various sin^ ^13 values from solar data (top) and 
solar+KamLAND data (bottom), before (lines) and after (colored regions) the SNO~salt data. 
The local minima in each plane after (before) SNO-salt data are marked by filled (open) dots. 



in the analysis of atmospheric and K2K data^, and to set Aml^-^i to infinity for the analysis 
of solar and KamLAND data. This implies furthermore that the effect of a possible Dirac 
CP- violating phase in the lepton mixing matrix can be neglected ^. We perform a general 
fit to the global data in the five-dimensional parameter space 5^2, sf^, Am^i, Am||, and 
show projections onto various one- or two-dimensional sub-spaces. 



B. Three flavour solar neutrino oscillations 

In this subsection we generalize the analysis of solar and KamLAND data presented in 
Sec. ini to three neutrino fiavours. Under the assumption of infinite Am^^^j only the new 
parameter 6'i3 appears for these data (see, e.g., Ref. |3J|). In Fig. HI we show the results 



of a three parameter fit (sin 6'sol, ArngQ^, sin ^^13) to solar and KamLAND data. Allowed 
regions are shown for various values of sin^6'i3 in the (sin^ 6'sol; Am^Q^) plane with respect 
to the global minimum, which occurs for sin^ ^13 = 0.02 including SNO-salt data, and for 



^ See Ref. l4J| for a two-mass_acale analysis of atmospheric data. 

^ The two Maiorana phases |43| do not show up in oscillations but do appear in lepton number violating 
processes 
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Figure 5: Ax^ profiles projected onto the sin'^ ^13 axis, for solar and KamLAND data, before 
(dashed lines) and after (solid lines) the SNO-salt experiment. 

sin^^i3 = 0.01 without the new data. Note that here we calculate the allowed regions at a 
given confidence level for 3 d.o.f.. From this figure one finds that the new SNO-salt data 
contributes to the disappearance of allowed regions when sin^ 613 increases. A comparison 
of the allowed regions before and after the SNO-salt data shows that their size is drastically 
reduced for all displayed values of sin^ 6*13. 

In Fig. we show the Ax^ projected onto the sin^ 613 axis for various data samples. 
From this figure one can observe an improvement on the constraint on sin^ ^13 thanks to the 
new SNO-salt data for values of sin^ 6^13 > 0.08. The shape of Ax^ can be understood from 
Fig. m From the upper panels of that figure one sees that increasing ^13 can be compensated 
to some extent by increasing ArngQ^- Since the new SNO-salt data disfavours large values of 
AmgQL the bound improves. Also KamLAND data acts in a similar way. Since the minimum 
around ArngQ^ ~ 7 x 10^^ eV^ is preferred the "jump" of the local minimum into the high- 
LMA-MSW region, which is visible in the lower panels of Fig. |3] and which leads to the 
kink in the Ax^ shown in Fig. |31 occurs at larger values of sin^^i3. We have also verified 
explicitly that assumptions regarding the statistical treatment of the solar boron fiux have a 
rather small effect on these results. Finally we note that, as will be seen in Subsec. IIIIDl the 
bound resulting from CHOOZ and atmospheric data still dominates the overall constraint 
on sin^ ^^13 in the global analysis. 
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parameter 


best fit 


2cr 


3cj 


5o- 


Am|i [lO-^eV^] 


6.9 


6.0-8.4 


5.4-9.5 


2.1-28 


Am|i [lO-^eV^] 


2.6 


1.8-3.3 


1.4-3.7 


0.77-4.8 


sin^ 6*12 


0.30 


0.25-0.36 


0.23-0.39 


0.17-0.48 


sin^ 6*23 


0.52 


0.36-0.67 


0.31-0.72 


0.22-0.81 


sin^ 6*13 


0.006 


< 0.035 


< 0.054 


< 0.11 



Table I: Best-fit values, 2cj, Scr and 5cr intervals (1 d.o.f.) for the three-flavour neutrino oscillation 
parameters from global data including solar, atmospheric, reactor (KamLAND and CHOOZ) and 
accelerator (K2K) experiments. 



C. Global analysis of solar, atmospheric, reactor and accelerator data 

In addition to the solar and KamLAND data described so far we now include all the 
charged- current atmospheric neutrino data, such as the most recent ones of the Super- 
Kamiokande [3, and MACRO experiments. The Super-Kamiokande data include 
the e-like and /z-like data samples of sub- and multi-GeV contained events (10 bins in zenith 
angle), as well as the stopping (5 angular bins) and through-going (10 angular bins) up- 
going muon data events. As previously, we do not use the information on Ur appearance, 
multi-ring fi and neutral-current events since an efficient Monte-Carlo simulation of these 
data would require a more detailed knowledge of the Super-Kamiokande experiment, and in 
particular of the way the neutral-current signal is extracted from the data. From MACRO 



we use the through 
analysis see Refs. 



i-g oing muon sample divided in 10 angular bins j27|. For details of our 
3lL |3J| and references therein. Furthermore, we include the first spectral 



data from the K2K long-baseline accelerator experiment js^. We use the 29 single-ring 
muon events, grouped into 6 energy bins, and we perform a spectral analysis similar to the 



one described in Ref. j46[. Fina^ 
the CHOOZ reactor experiment 



■ly, we take into account in our fit the constraints from 
29|. The best fit values and the allowed intervals of the 
parameters 5^2, -5235 -Sis, Am^i, Am|^ from the global data are given in Tab. HI This table 
summarizes the current status of the three fiavour neutrino oscillation parameters. 

For completeness we show in Fig. IHlthe projection of the allowed regions from the global 
fit onto the plane of the atmospheric neutrino parameters. We find that the first 29 events 
from K2K included here start already to constrain the upper region of Am^^^. We note 



that recently the Super-Kamiokande collaboration has presented a preliminary reanalysis of 
their atmospheric data ^. The up-date includes changes in the detector simulation, data 

analysis and input atmospheric neutrino fiuxes. These changes lead to a slight downward 
shift of the mass-splitting to the best fit value of A?ti^ = 2 x 10~^ eV^. Currently it is 
not possible to recover enough information from Ref. j43] to incorporate the corresponding 
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Figure 6: Projections of the allowed regions at 90%, 95%, 99%, and 3a C.L. for 2 d.o.f. onto the 
plane of sin^ 6'atm = sin^ 623 and Am^^^j = Am|^ . The regions delimited by the lines correspond 
to atmospheric+solar+CHOOZ data, for the colored regions also K2K and KamLAND data is 
added. Also shown is the Ax^ as a function of sin^ ^atm and Am^^^^, minimized with respect to all 
un-displayed parameters. 

changes in our codes. We note, however, that the quoted value for Am^^^ is statistically 
compatible with our result. For Am^^^ = 2 x 10~^ eV^ and maximal mixing we obtain a 
Ax' = 1.3. 



D. Status of and the mass hierarchy parameter a = Amg^^j/Am: 



atm 



In this subsection we discuss in detail the current information on the small parameters 
relevant for sub-leading oscillation effects. First we consider the mixing angle ^13, which at 
the moment is the last unknown angle in the three-neutrino mixing matrix. Only an upper 



2c 



Currently 



38|,|48|,|49|). 



bound exists, which used to be dominated by the CHOOZ reactor experiment 
a large effort is put to determine this angle in future experiments (see, e.g., Refs. 

In Fig. [7|we show the Ax^ as a function of sin^ ^13 for different data sample choices. One 
can see how the bound on sin^ Oi^ as implied by the CHOOZ experiment in combination 
with the atmospheric neutrino data still provides the main restriction on sin^^i3. We find 
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Figure 7: profiles projected onto the sin^ ^13 axis, for solar +KamL AND, atmo- 

spheric+CHOOZ, and for the global data before (dashed lines) and after (solid lines) the SNO-salt 
experiment. 

the following bounds at 90% C.L. (3cr) for 1 d.o.f.: 

{ 0.070 (0.12) (solar+KamLAND) 
0.028 (0.066) (CHOOZ+atmospheric) (5) 
0.029 (0.054) (global data) 

However, we note that the solar data contributes in an important way to the constraint 

on sin^ £12. for lower values of Am^^^. In particular, the down- ward shift of Am^^^,, reported 

,2 , 



in Ref. |42[ implies a significant loosening of the CHOOZ bound on sin 613, since this bound 
gets quickly weak when Am^^^^ decreases (see, e.g., Ref. 5^). Such loosening in sensitivity 
is prevented to some extent by solar neutrino data. In Fig. |H1 we show the allowed regions 
in the (sin^ ^13, Am^^^J plane from an analysis including solar and reactor neutrino data 
(CHOOZ and KamLAND). One finds that, although for larger Am^^^„ values the bound on 
sin^ 6'i3 is dominated by the Chooz + atmospheric data, for low Am^^j^j = 10~^eV^ the solar 
+ KamLAND bound is comparable to that coming from Chooz + atmospheric. 

For example, fixing Am^^^ = 2 x 10^'^ eV^ we obtain at 90% C.L. (3a) for 1 d.o.f. the 
bound 

sin^ ^13 < 0.035 (0.066) (Am^^^, = 2 x 10"^ eV^) , (6) 
which is only marginally worse than the bound from global data shown in Eq. 
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Figure 8: Allowed regions in the (sin^ 613, Am^^M) plane at 90%, 95%, 99%, and 3a from CHOOZ 
data alone (lines) and CHOOZ+solar+KamLAND data (colored regions). 

For the exploration of genuine three-flavour effects such as CP-violation the mass hi- 
erarchy parameter a = Am'^^^/ Aml^-^^ is of crucial importance since, in a three-neutrino 
scheme, CP violation disappears in the limit where two neutrinos become degenerate jsil ]. 
Therefore we show in Fig. IHl the Ax^ from the global data as a function of this parameter. 
Also shown in this figure is the Ax^ as a function of the parameter combination asin26'i2, 
since to leading order in the long baseline Ue z/^ oscillation probability solar parameters 
appear in this particular combination We obtain the following best fit values and 3a 
intervals: 

a = 0.026, 0.018 < a < 0.053, 

- - ' (7) 
a sin 2^12 = 0.024 , 0.016 < a sin 2^12 < 0.049 . 



IV. CONCLUSIONS 



We have performed a global analysis of neutrino oscillation data in the three-neutrino 
scheme, including the recent improved measurement of the neutral current events at SNO. 
We include in our fit all current solar neutrino data, reactor neutrino data from KamLAND 
and CHOOZ, the atmospheric neutrino data from Super-Kamiokande and MACRO, as well 
as the first data from the K2K long-baseline accelerator experiment. We have discussed 
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the implications of the improved neutral current measurement on the determination of the 
solar neutrino oscillation parameters. The SNO-salt results reject the previously allowed 
high-mass branch of AmgQ^ cit about 3a. Moreover, it rules out maximal solar mixing at 
more than 5a, precluding the possibility of bi-maximal models of neutrino mixing. 

Furthermore, we have determined the current best fit values and allowed ranges for the 
three-flavour oscillation parameters from the global oscillation data. Our three-neutrino re- 
sults are summarized in Tab.|l]and in Fig.^J showing the allowed regions and projections 
for the five oscillation parameters sin^ 6'sol5 sin^ ^'atm, sin^ 6*13, ArrigQ^, Am^^^^. In addition we 
have discussed in detail the limits on the small mixing angle sin^^is, as well as the hierar- 
chy parameter a = ArngQ^/^^ATM- These small parameters are relevant for genuine three 
flavour effects, and restrict the magnitude of leptonic CP violation that one may potentially 
probe at future experiments like super beams or neutrino factories. In particular, we have 
seen how the improvement on the sin^ ^13 limit that follows from the new SNO-salt data 
solar neutrino experiments can not yet match the sensitivity on sin^ 613 reached at reactor 
experiments. However, the solar data do play am important role in stabilizing the constraint 
on sin^ ^^13 with respect to variation of Am^^^j. For small enough Aml^^ values the solar 
data probe sin^ 6*13 at a level comparable to that of the current reactor experiments. 
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